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Abstract 

This paper introduces in full detail a methodology 
for the measurement of neutron yields and the nec- 
essary efficiency calibration, to be applied to the in- 
tensity measurement of neutron bursts where indi- 
vidual neutrons are not resolved in time, for any 
given moderated neutron proportional counter array. 
The method allows efficiency calibration employing 
the detection neutrons arising from an isotopic neu- 
tron source. Full statistical study of the procedure 
is descripted, taking into account contributions aris- 
ing from counting statistics, piling-up statistics of real 
detector pulse-height spectra and background fluc- 
tuations. The useful information is extracted from 
the net waveform area of the signal arising from the 
electric charge accumulated inside the detector tube. 
Improvement of detection limit is gained, therefore 
this detection system can be used in detection of low 
emission neutron pulsed sources with pulses of dura- 
tion from nanoseconds to up. The application of the 
methodology to detection systems to be applied for 
D — D fusion neutrons from plasma focus devices is 
described. The present work is also of interest to the 
nuclear community working on fusion by magnetic 
confinement and lasers, and working on neutron pro- 
duction by accelerators or similar devices. 
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1 Introduction 

Foil activation is usually used as a standard tech- 
nique for neutron yield measurements in pulsed fu- 
sion sources. From the early plasma focus (PF) re- 
search [1, 2], activation of Silver and Indium foils 
have been used [3, 4, 5], although the so called silver 
activated Geiger counter is the most known detector 
from both. In PF devices [6], depending on the filling 
gas, neutrons from D — D reactions or D — T reac- 
tions are produced with typical energies of 2AbMeV 
and li.lMeV respectively. For fast neutrons, acti- 
vation detectors usually require neutron moderation 
to thermal energies. Notwithstanding, when neutron 
intensity is high enough, typically higher than 10®n 
per source pulse, activation by fast neutrons could 
be used as in the case of the Indium or Beryllium 
counters. [5, 7, 8]. Detectors based on activation 
by thermal neutrons have usually detection limits 
higher than lO^n per source pulse [3, 5, 9]. 

In recent years, several groups around the world 
have started research programs on small scale low 
energy plasma focus devices, especially in the sub- 
kiloJoule energy range [10, 11, 12, 13, 14, 15, 16, 
17, 18, 19] and in the sub-hundredJoule energy 



range [20, 6, 21, 22, 23, 24]. As a consequence, re- 
quirements on detector sensitivity were increased in 
order to characterize such devices. For the lov^r emis- 
sion regime (Y < 10'' per source pulse) it is possible 
to use moderated gas proportional counters based 
on ^He or BF3 filled tubes [25]. In these detectors 
thermal neutron capture occurs with cross sections 
one order of magnitude higher than Indium or Silver 
activation cross section. Therefore, detection limits 
in neutron proportional counters should be at least 
one order of magnitude lower The disadvantage of 
these detectors for pulsed sources is the pulse piling- 
up produced during high instantaneous count rates 
due to the burst of fast neutrons, which make it im- 
possible to count single neutron pulses by standard 
nuclear electronics. To deal with this problem two 
schemes of use have been proposed: 

i) Use of high counting rate electronics and con- 
trolling of the instantaneous count rate by 
detector-source separation. This technique has 
been reported to be successful in resolving yields 
higher than lO^n per source burst [26]. 

ii) Cross calibration on the accumulated charge 
(preamplifier output signal) generated in the 
counter tube by the neutron burst using as neu- 
tron reference a foil activation detector [27, 16] . 
The comparison of the response of a reference 
detector with the proportional counter allows 
to obtain a calibration factor in the high emis- 
sion regime. Thus, by means of extrapolation, 
it is possible to detect yields of the order of 
lO'^ — lO'^n/burst in the low emission regime us- 
ing the moderated proportional counters. Due 
to the fact that detection limits are transferred 
by the calibration process from the reference de- 
tector, the uncertainties in single measurements 
are of the order of the yield or greater Thus, 
no improvement in measurement uncertainties 
is obtained for the low emission regime despite 
the high efficiency of neutron proportional coun- 
ters tubes in comparison with foil activation de- 
tectors. However, for bursts of the order of 
lO'^ij/burst, it is still possible to measure the 
pulsed neutron emission rate with an accuracy 



of the order of 40% by taking averages on large 
samples of shots [21]. 

The current state of the art in the use of propor- 
tional counters for pulsed fast neutrons lacks a cal- 
ibration technique to allow reducing detection lim- 
its below lO^n per source burst in single shot mea- 
surements. To overcome this deficiency, in the fol- 
lowing, a methodology is proposed for the analysis 
of the detector output signal to find the number of 
piled-up individual detected events per burst and its 
uncertainties, while allowing for detector efficiency 
to be established through the measurement of neu- 
trons arising from an isotopic neutron source. Thus, 
it is possible to obtain the real neutron yield. The 
detailed study of the pulse piling-up statistics and 
other sources that influence the measurement pro- 
cess, namely counting statistics and background, are 
described in what follows. 

2 Working principles for fast 
neutron counting through 
a moderated proportional 
counter tube 

In the so called continuous regime, as depicted in fig. 
1, detection system consists of a moderated counter 
tube connected to a charge preamplifier and an am- 
plifier, and finally connected to a counting system. 
As a consequence from the detection process in the 
counter tube, charge is produced and collected at the 
input terminal of the preamplifier The output of the 
preamplifier is a signal whose pulse height is propor- 
tional to the integrated charge provided at the input. 
Decay of the preamplifier output is characterized by 
an exponential fall and RC time constant. Once in 
the amplifier, typically a linear amplifier, the signal 
is shaped and amplified in order to make it suitable 
for counting systems. To allow proper counting, the 
count rate should in general satisfy ?- << (RC)^^. 
In fig. 1 it is also shown the output in the so called 
pulsed regime. In this case the preamplifier output 
is directly connected to a digital oscilloscope. When 
the detection system is irradiated by an extremely 
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Figure 1: Detection system setup and signal forma- 
tion stages for (a) Continuous regime and (b) Pulsed 
regime. 



short pulse of fast neutrons (r < lO^ns), the pulse is 
spread over tens of microseconds as a result of neu- 
tron moderation to thermal energies. When fluence 
is high enough to make the instantaneous count rate 
to be Tins » (RC)^^, single neutron signals will be 
piled-up thus smoothing the output signal and mak- 
ing it impossible to count individual detected events 
by standard nuclear electronics. In contrast for lower 
fluence, when rms < {RC)~^, signal shape is not 
smooth and eventually it is possible to distinguish 
single neutron pulses from the output signal. In both 
cases notwithstanding, the total number of detected 
events should be proportional to the output signal 
area. Even in the case of pulse piling-up, the num- 
ber of detected events could be always obtained only 
from the net signal area, provided that space charge 
accumulation in the detector tube does not impair 
linearity. 

When detection system is set up for pulsed regime 
{"charge integration mode"), there are three sources 
of fluctuations that influence the measurement pro- 
cess: 

i) Counting statistics: Fluctuations arise from the 
detection process and are well described by the 
Poisson distribution. 

ii) Pulse piling-up statistics: Arise from the ''wall 
effect" taking place inside the counter tube, 
which results in a very asymmetrical distribu- 
tion. Therefore, when pulse piling-up happens, 
the sum of two or more single random detected 
events generate overlapping distributions which 
give rise to differences between the true and the 
expected number of detected events. 



iii) Background: The extracted signal from the de- 
tection system corresponds to the total charge 
generated in the tube, which is equivalent to the 
integral of the waveform acquired by the oscillo- 
scope. Thus, electrical noise acts as background. 
Besides electrical noise, natural background is 
always present although in some cases it could 
be negligible. 
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3 Detection system and experi- 
mental setup 

Detection systems available in our facilities were 
reported previously by Moreno et al [27] . They are 
based on a paraffin wax moderator (45 x 15 x 15cm^) 
and a ^He tube (model LND 2523) sheathed by a 
3mm thick lead sheet to stop x-rays and prevent their 
subsequent detection. The moderator is surrounded 
by an external cadmium sheet to absorb thermal 
neutrons. The tube is connected to a preamplifier 
(CANBERRA 2006) and the high voltage feed is 
set to 4fcy. In order improve the signal-to-noise 
ratio, the preamplifier is adjusted with a conversion 
gain of 235mV/M — ion — par which produces a 
5x voltage output factor scale [28]. The RC time 
constant is set to 50/^is, this is the fall time for 
individual pulses when the input pulse risetime is 
less than hundred nanoseconds. When used for the 
pulsed regime, signal output from the preamplifier is 
recorded by a digital oscilloscope (Tektronix model 
TDS 684) at IMfl input impedance, AC coupling, 
and 20MHz bandwidth. Horizontal scale is set in 
accordance with moderator decay time (governing 
neutron feeding into the "^He tube), in the case of 
our detectors time scale is set to 200fis/div. Direct 
measurements of the total signal area are obtained 
using the internal integration oscilloscope function 
on a 1500/is time window, which starts close to the 
trigger point. The oscilloscope is triggered by the 
electromagnetic pulse generated from the electrical 
discharge. 

Detector efficiency or calibration factor (Jc) for fast 
neutrons was obtained operating in the continuous 
regime by comparing neutron intensity from a stan- 
dard isotopic source, used as neutron reference, and 
the registered net counting rate. The last is valid in 
the approximation that after many neutron scatters 
in the moderator the result is highly independent 
of the neutron input spectra; that is, the standard 
source spectra or that originated at the Plasma Focus. 
To satisfy this approximation the fast neutron refer- 
ence should have a neutron energy spectra as similar 
to the fusion spectra as possible. Although D — D 



reaction at low incident energies produces 2.45Mey 
neutrons, observed fusion spectra around 2AbMeV 
and extending from IMeV to 3.5MeV have been 
reported in the literature [29, 5]. To fulfill this re- 
quirement a certified ^''^Cf isotopic source was used. 
This source is characterized by a mean energy be- 
tween 2.1 — 2.3MeV and its spectrum extends from 
0.1 — lOMeV. Source neutron emision at the time 
of experiments was 6.5 x lO^s"^. Measurements 
of count rate were carried out at different positions 
from 14cm to 170cm. For this purpose, the preampli- 
fier output was connected to a Tennelec TC-244 am- 
plifier and then to a multichannel analyzer Camberra 
Multiport II controlled by the interface GENIE2K. The 
gross counting was corrected by discounting natural 
background and air attenuation according to that re- 
ported by Eisenhauer et al. [30]. To account for 
solid angle, the neutron intensity was corrected by 
the source anisotropy factors [31], which were pre- 
viously characterized [32], and thus the number of 
neutrons arriving to the detector front face was cal- 
culated using the formula provided by Gotoh et al. 
for a rectangular slit [33, 34]. These same measure- 
ments of count rate provided the pulse height spec- 
trum of each proportional counter. In the following 
sections results are referred to one of our detection 
systems which is called ■^He-206 system. 

4 Piling-up statistics for propor- 
tional counters 

Piling-up statistics was studied by means of the 
acceptance-rejection Monte Carlo method, which 
was used to numerically construct pulse piling-up 
distributions. The starting point was to obtain the 
waveform area spectrum for single events. Due to 
pulse height at the preamplifier output being propor- 
tional to the accumulated charge provided at the in- 
put, then for single events waveform area is given 
hy Ap = Vp ■ FA, where Vp is the pulse height and 
FA a constant. Thus, the pulse area spectrum could 
be obtained from the pulse height spectrum through 
rescaling it by the constant FA. FA depends only on 
the detection system. In order to characterize FA, 
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Figure 2: Diagram of a typical output signal from 
the detection system associated to a single detected 
event. 



Table 1: Characteristic preamplifier parameters for 
single neutron events acording to eq. 1 in •^He-206 
detection system. 



Parameter 




Value 






ri(/xs) 


1.24 


± 


0.01 


(0.5%) 


T2(/XS) 


75.1 


± 


0.3 


(0.4%) 




17.7 


± 


0.1 


(0.6%) 


FAiiis) 


91.6 


± 


0.3 


(0.3%) 



let us first consider a typical output signal from the 
preamp associated to detection of a single event, as 
shown in fig. 2. This signal can be modeled by the 
mathematical expression 

(b ,ift<-toi. 
Sit) = I Vpil - e-(*+*«i)M) +b , if - toi < i < to2- 

(1) 

Thus the signal area is given by 

Ap = Vp- (t2 + ri(e-(*°^+*°=)/^i - 1) + toi + ^02) 
where 

i^A = (T2+Ti(e-(*°i+*°^)/^i -l)+toi+to2). (2) 

To obtain proper values to calculate FA, hundreds 
of waveforms from single neutron events were an- 
alyzed using the model of eq. 1 and the non-linear 
least squares (Levenberg-Marquardt method) routine 
provided by the function fit of the GNUPLOT pack- 
age [35]. Results of this analysis are presented in 
table 1. Using these results it was possible to con- 
struct the probability density function (pdf) and the 
cumulative distribution function for the proportional 
counter tube as shown in figure 3. 

Piling-up statistics was studied by simulating the 
random detector response in terms of the waveform 




Signal area (v ■ ms) 

Figure 3: Probability density function (pdf) and cu- 
mulative distribution function for the proportional 
counter tube in •^He-206 detection system. 
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area and the other statistics associated to the mea- 
surement process. Let m be the total number of de- 
tected events and x() a random number extracted 
from the counter tube pdf (fig. 3) by means of the 
acceptance-rejection Monte Carlo method. When m 
events are detected, the random ideal response of the 
detection system is given by 



Xrim) 



xi). 



(3) 



Thus the probability density function from the pulse 
piling-up statistics is obtained as a function of m by 
sampling of equation 3. 

Samples were generated numerically for m : 
1, . . . , 1000. The sample size choice was 20000. Re- 
sults from the '^He-206 detection systems are shown 
in fig. 4. As expected, for m ~ 1 the numerical algo- 
rithm reproduces the experimental pdf of the counter 
tube. From numerical results it is observed that when 
m < 50 pulse piling-up distribution is highly asym- 
metrical, for ?n > 100 asymmetry is still present at 
tails of the distributions, finally when rn > 500 tail 
differences are not apparent from histograms. From 
the central limit theorem, the pulse piling-up statistic 
is expected to converge to the normal distribution. 
This property was evaluated at the 1% significance 
level by means of Kolmogorov-Smirnov and Jarque- 
Bera normality test. Convergence of the counter tube 
pdf (fig. 3) to the normal distribution was found 
in terms of cumulative probability (K-S test) when 
TO > 90, while in terms of asymmetry and peaked- 
ness (J-B test) when to > 800. 



5 Background readings 

False positive readings are generated in the detection 
system by two mechanisms: 

i) Natural background: Corresponds to events 
from the environment which are characterized 
by a constant count rate. Natural background 
count rate in our experimental room is usually 
found in the range 0.3 — 1.0s~^. These measure- 
ments were done under the continuous regime. 




0.000 
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Figure 4: Pulse piling-up distributions from sampling 
of eq. 3 for ''He-206 detection system. 



For the pulsed regime the measurement time 
window is 1500/is, therefore contribution to the 
total number of detected events by natural back- 
ground is negligible. 

ii) Electrical background noise (EBN) : In absence 
of detected events, the electrical noise con- 
tributes randomly to the total waveform area 
(detector reading), thus it behaves as a back- 
ground reading. When the operative parameters 
in the detection system are fixed, readings from 
EBN respond to a normal distribution. The lat- 
ter is concluded from analysis of EBN samples 
by the Kolmogorov-Smirnov test at 10% signifi- 
cance level. 

These facts imply that the gross detector read- 
ing (XTgross) should be corrected to discount back- 
ground contributions. Let B and dB respectively be 
the sample mean and sample standard deviation as- 
sociated to the EBN. Therefore, the net detector read- 
ing is obtained from 

^Tnet = ^Tgross - B. (4) 

For a specific time window, the EBN sample sta- 
tistical parameters depend on oscilloscope horizontal 
and vertical scale, they are affected by changes in lab- 
oratory room temperature and should be monitored 
in-situ together with neutron measurements for shots 
with no neutron detection. 

6 Counting model from the 
waveform area 

To arrive at a counting model as a function of the 
net waveform area, let us first consider a sequence 
of random values associated to m' detected events: 
{Xti,Xt2, ■ ■ ■ , Xtu}- Therefore the following ex- 
pressions are satisfied 

m' = uiXti 
m! = 02X7^2 

m' = anXrn- 



Table 2: Linear fit results using model from equation 
6. 



Parameter Value 

Range m : 1 - 300 

Sample size 20000 

6771 1 

a{V-^ ■ ^ls-^) 0.44459 ± 0.00001 
fit r2 0.99999907 



Selecting = a as a constant, let the following ?ti' 
be defined as follows, 

7n'i = a ■ Xti 

771^ = a ■ Xt2 

= a-Xxn- 

By properly choosing a, convergence of the distribu- 
tions ensures that 

E{77i'^) = E{a ■ Xt.,) = a ■ E{Xti) = m'. (5) 

Being E{x) the mathematical expectation of the vari- 
able X. This constant a was evaluated using samples 
from the pulse piling-up statistic. Detected number 
of events was plotted against the sample mean of the 
waveform area, 

777{E{XTi)) = a-E{XT^). (6) 

Thus, a linear least squares through the origin yields 
the value of a. Results of this analysis are shown in 
table 2. 

Finally the counting model is defined by means of 
the rounding function, which round the argument to 
the nearest integer, as follows 

m„,, = G{XT)^[a-XTl (7) 

where m,„c corresponds to a predicted value for the 
number of detected events from the signal wave- 
form area or accumulated charge. Once the count- 
ing model has been obtained, it becomes necessary 
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to study the predictive power of the method and the 
uncertainties. As aheady stated in section 2, in prac- 
tical situations fluctuations are due to the counting 
statistics, which responds to the Poisson distribution, 
also due to the pulse piling-up statistics whose dis- 
tribution has been studied in section 4 and to the 
electrical background noise which is governed by the 
normal distribution (see section 5) . 

Predictive power was studied by simulating ran- 
dom measurements as predictions from the count- 
ing model. For that purpose pseudo-random num- 
bers from the Poisson and normal distributions 
where required. Due to its proven mathematical 
properties although being slower than others, the 
gsl_rng_ranlux389 pseudo-random number gener- 
ator provided by the GNU Scientific Library (GSL) 
was used [36] . The algorithm for simulating random 
measurements and model predictions is described as 
follows: 

1. Let m be the true number of detected events. 

2. Let be a pseudo-random integer extracted 
from the Poisson distribution with expected 
value 771. 

3. Let Xtec be a pseudo-random number extracted 
from the pulse piling-up statistic associated to 
rriec detected events. 

4. Let Bran be a pseudo-random number extracted 
from a Gaussian distribution Qi ~ 0) with a = 
dB. Thus a new variable is defined by Xxnet ^ 
Xxec + Bran- At this poiut, Xxnet corrcspouds 
to the simulated net detector reading. 

5. Finally a prediction from the counting model is 
obtained by m,„c = G{XTec)- 

Sampling of the above described algorithm al- 
lows further calculation of the statistical parame- 
ters Eim^c), E{777,nc), E{Bran), ^(mec), (j{m,nc) and 
a{Bran), ^nd comparison of them with the true num- 
ber of detected events. 

Samples were created for the interval 1 < m < 
300 with a step size of 5m = 2. Sample size choice 
was 5000. For each value of m, dB was varied in the 
interval < dB < IQ V ■ us with steps of 5dB = 



0.25V ■ IIS when dB < 2.0V ■ ^s, and SdB = 0.5V ■ ^s 
when 2.0 < dB < lOV ■ /is. 

Analysis of numerical results to determine predic- 
tive power and model uncertainties was done by 
means of linear models and hypothesis testing. For 
predictive power, if predictions from the counting 
model iE{m,nc)) converge to the true number of de- 
tected events (to), then these two variables should 
be related in a X — Y plot by the identity function. 
To evaluate shift of the counting model caused by its 
definition or by numerical reasons, a general linear 
model is used: y{x) = Pix + Pq. Thus, the hj^othesis 
to be tested by the F-test is given as follows: 

iJo : /3i = 1 y /3o = Null hypothesis 



Ha:f3iy^l 



versus 

y 



Alternative 
hypothesis 



If null hypothesis is accepted for a given signifi- 
cance level a, then a high predictive power is con- 
cluded for the counting model. Otherwise, the count- 
ing model is assumed to be poorly predictive and 
should be rejected. 

For the F-test, the significance level was chosen to 
be 1% (a = 0.01). Therefore, our sample with 148 
degrees of freedom has a F-test critical value of 4.751. 
At the top of figure 5 results are shown of applica- 
tion of the F-test on samples as a function of dB for 
the study of counting model predictive power. For 
each sample, it was verified that the expected value 
of the counting statistics also converge to the true 
number of detected events. From figure 5, it is ob- 
served that all calculations are less than the critical 
F-value, thus the counting model defined by equation 
7 is concluded to be highly predictive, this is to say 



E{mmc) = m. 



(8) 



The last being valid at least in the interval under 
study: to < 300 and dB < lOV ■ ^s. 

Once counting model is accepted as predictive, it is 
necessary to study uncertainties of the model. With 
this purpose, two cases were considered: 
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Figure 5: Test of hypothesis results for predictive 
power of counting model and uncertainties. Criti- 
cal value for the F-test at the 1% significance level is 
shown as a solid line. 

• Case 1: Uncertainties of the counting model 
due only to counting statistic and pulse piling- 
up statistic. This case is equivalent to consider 
sample generation when 

lim Bran 0. (9) 

• Case 2: Uncertainties of the counting model due 
to counting statistic, pulse piling-up statistic and 
electrical background noise. 

For case 1 samples in the interval 1 < m < 300, 
with 6m. ~ 1 and sample size of 20000, were ana- 
lyzed. As a result, it was found that variance of pre- 
dictions from the counting model are proportional to 
the expected value for predictions from the counting 
model, 

o-'^{mrnc) = /o • E{m,nc) = fo ' m. (10) 

Thus, distributions generated by the counting 
model, neglecting background effects, could be said 
to be Poisson-like because effects of the counting 



statistic dominate over pulse piling-up statistic. Last 
conclusions are supported by numerical results as 
shown in figure 6 and table 3. On one hand, samples 
from the counting statistics were analyzed by linear 
least squares (LLS) using the model a{mec)^ = ^ -m. 
Use of the T-test of hypothesis allowed to conclude 
for the 1% significance level that 7 = 1, as expected. 
On the other hand, analysis of sample data by LLS 
and model results from equation 10 allow to con- 
clude that /o = 1.065. Therefore, it is possible to say 
that pulse piling-up statistical contribution is 6.5% 
the contribution from the counting statistic to count- 
ing model variance when neglecting background ef- 
fects. 




m, Number of detected events 

Figure 6: Standard deviation modeling from sample 
data for ■^He-206 detection system. 



Table 3: Counting model uncertainties study for case 
1. 





Parameter 


Value 


Counting 


fo 


1.065 ± 0.0006 


model 


2 

7' 


0.9995687 


Counting 


7 


0.9998 ± 0.0006 


statistic 


r2 


0.99958056 



For case 2 is postulated that counting model vari- 
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ance is given by, 

cr^ini„ic) = fa ■ E{mec) + cr^{a ■ Bran) 



Table 4: Characteristic values for jc and /47r in the 
interval 14cm < r < 40cm. 



fo-m + a^ ■ dB'^ 



(11) 



sample ■ 



To demonstrate the last hypothesis, predictive 
power of equation 1 1 was studied. If the right side 
of equation 11 converges to the left one, then the 



plot (7^ {rrimc) v/s fo ■ni + ■ dB 



sample 



should be 



given by the identity function. By considering sam- 
ples with the same characteristics than those set up 
for the study of counting model predictive power, it 
was obtained that equation 1 1 is highly predictive for 
the sample counting model standard deviation, as it 
could be verified at bottom of figure 5. It is impor- 
tant to note that to achieve predictability, the sample 
standard deviation of the background electrical noise 
should be used. 

Summarizing results so far, it has been demon- 
strated that definition of the counting model by equa- 
tion 7 ensures E{mmc) — "m and u^{m„ic) = /o • + 

' ^^^ample ■ 

7 Neutron yield measurements 
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Figure 7: Mean relative error for ^He-206 detection 
system. 

Let us consider a single shot from the pulsed 
source, e.g a plasma focus device. Now we deal with 
the problem of estimating the expected number of 
detected events and its uncertainties. Let Xxnet the 



Parameter 


Mean 




14.2 ±0.1 


/4x = 4:TT/n 


20.45 



net detector response (see equation 4). Then, for 
this single measurement the best estimation of the 
expected number of detected events comes from the 
measurement itself [37], 



G{XTnet) = [a ■ Xxnet]- (12) 



Thus, the uncertainty associated to the expected 
number of detected events in a single measurement 
is calculated from 



fo-m + a? ■ dB^ 



sample 



= Jf0-GiXTnet)+Ci^-dBl 



sample 



(13) 



/o • [a ■ XTnet] + ■ dBl 



sample 



Once obtained an estimation of the number of de- 
tected events has been obtained, the neutron yield is 
easily deduced by using the calibration factor of the 
detection system (j^). Consequently, the number of 
neutrons per solid angle is given by 



Y =jc • [a ■ Xtu 




Let O be the solid angle subtended by the detector 
at its measurement position, then f^j, = An /ft. Thus 
under the assumption of isotropic emission, the total 
neutron emission per pulse is given by 



YiTT — fin ■ Y- 



(15) 



In order to study detection limits, let us consider 
table 4 where characteristic values for jc and 
of the typical positioning interval in an extremely 



10 



low energy plasma focus device (PF-50J, [21]) are 
shown. In addition, from [a ■ Xxnet] = ^/jc and 
Y = Ki7r//47r the standard relative error is then cal- 
culated by 

%dY(Y) = 100- 



/o • jc • fi-rr f ■ dBsample ' jc ' /47r 




47r 



(16) 



Plots from equation 16 are shown as a function 
of dBsample in figurc 7. It is observed from this fig- 
ure that when total neutron emission is higher than 
IQ^n/pulse, uncertainties are less than 10%. For 
IQ^n/pulse uncertainties are of the order 20%. Fi- 
nally, it has to be pointed out that reasonable un- 
certainties (< 40%) are still obtained in the range 
2 • 10^ < 14^ < IQ^n/ pulse. The last is a direct con- 
sequence of the methodology proposed in this work. 



8 Concluding remarks 

The statistics of the signal generated from pulse 
piling-up (pulse piling-up statistics) in a neutron pro- 
portional counter when used in the herein called 
pulsed regime has been studied. It has been found 
that pulse piling-up statistics converge to a nor- 
mal distribution typically when the number of piled- 
up pulses is higher than a hundred. From the 
pulsed piling-up statistics, the counting statistical be- 
haviour and background contributions, a measure- 
ment methodology for bursts of neutrons has been 
developed. The method is based on a counting model 
for the number of detected events based upon the 
net accumulated charge at the output of the detec- 
tion system (waveform area). It has been demon- 
strated that the model, as defined in equations 7 and 
10, ensures E{mmc) = m and C7^{mmc) — fo ■ + 



■dBi 



, when m is the true number of detected 



sample' 

events. Then by using equation 14 it is possible to 
calculate the neutron emission, bearing in mind that 
the model holds while space charge accumulation in 
the detector tube does not run so high as to impair 
linearity. By this methodology, it has been possible 



to achieve detection limits almost two orders of mag- 
nitude lower than those from state of the art tech- 
niques. The methodology is quite general and could 
be reproduced by other groups in order to calibrate 
neutron detectors based on moderated proportional 
counters for fast neutrons. 
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